function res=Anz(lambda_t,x_t,n,xi,chi,mu_l,alpha,v,mu_x,rho,phi,psi,omega,etaz,piz,ngridpoints)

b=fliplr(piz*psi*(1-omega.^(0:n-1))/(1-omega)+etaz);

res=compute_exp(lambda_t,x_t,n,b,xi,chi,mu_l,alpha,v,mu_x,rho,phi,ngridpoints);